Transport and elastic scattering times as probes of the nature of impurity scattering 

in single and bilayer graphene 
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Both transport rtr and elastic r e scattering times are experimentally determined from the carrier 
density dependence of the magnetoconductance of monolayer and bilayer graphene. Both times and 
their dependences on carrier density are found to be very different in the monolayer and the bilayer. 
However, their ratio Tt r /r e is found to be close to 1.8 in both systems and nearly independent of the 
carrier density. These measurements give insight on the nature (neutral or charged) and range of the 
scatterers. Comparison with theoretical predictions suggests that the main scattering mechanism 
in our samples is due to strong (resonant) scatterers of a range shorter than the Fermi wavelength; 
likely candidates being vacancies, voids, ad-atoms or short-range ripples. 

PACS numbers: 63.22.Np, 73.23.Hk, 73.21.Ac 



Since the discovery of the fascinating electronic prop- 
erties of graphene pQ due to its electronic spectrum with 
linear dispersion and a perfect electron-hole symmetry at 
the Fermi level [2], the nature of defects has been shown 
to play an essential role in determining the carrier den- 
sity (n c ) dependence of the conductance. The wavevec- 
tor and energy dependences of the impurity potential are 
known to determine the characteristic scattering times of 
the carriers. It is important to distinguish the transport 
time Tt r , which governs the current relaxation and enters 
the Drude conductivity (a), from the elastic scattering 
time r e , which is the lifetime of a plane wave state [3]. 
Since r tr and r e involve different angular integrals of the 
differential cross section, they differ as soon as the Fourier 
components of the potential depend on the wavevector q. 
A large ratio r^ r /r e indicates that scattering is predomi- 
nantly in the forward direction, so that transport is not 
affected much by this type of scattering. This is the case 
in 2D electron gases (2DEG) confined to GaAs/GaAlAs 
hetero junctions with the scattering potential produced 
by remote charged Si donors [4], where r tr /T e is found to 
be larger than 10. 

The nature of the main scattering mechanism limit- 
ing the carrier mobility in graphene is still subject to 
controversy. It has indeed been shown [5HZ] that "white 
noise" (q independent) scattering leads to a weak (loga- 
rithmic) dependence of <r(n c ), in contradiction with ex- 
periments which typically find a linear increase. In con- 
trast, scattering on charged impurities originates from 
a q dependent screened Coulomb potential described in 
the Thomas Fermi approximation [8HTD]. This leads to 
a linear cr(n c ) both for a monolayer (ML) and a moder- 
ately doped bilayer (BL). Recent experiments performed 
to probe this question measured the change in a upon 
immersion of graphene samples in high-i\~ dielectric me- 
dia. Their conclusions differ [IT] . Alternate explanations 
involve resonant scattering centers with a large energy 



mismatch with the Fermi energy of carriers [7J [12] . 

In order to gain insight into the scattering mechanism 
in graphene, we have extracted r e and r tr from mag- 
netotransport in monolayer and bilayer graphene sam- 
ples. In high magnetic field, when the cyclotron fre- 
quency is larger than 1 /r e , the magneto-conductivity ex- 
hibits Shubnikov de Haas (ShdH) oscillations related to 
the formation of Landau levels. The broadening of these 
levels at low temperature yields r e , while the low field 
quadratic magneto-conductivity yields r tr • 

The samples were fabricated by exfoliation of natural 
graphite flakes and deposition on a doped silicon sub- 
strate with a 285 nm thick oxide. The carrier density 
can be tuned from electrons to holes through the charge 
neutrality point by applying a voltage on the backgate. 
The ML and BL samples were identified using Raman 
spectroscopy. The electrodes were fabricated by electron 
beam lithography and either sputter deposition of 40 nm 
thick palladium (samples A and B), or Joule evaporation 
of a bilayer 5nmTi/70nm Au (other samples C,D and E). 
We mostly discuss samples A and B, a ML and a BL 
of respective dimensions W = 1.6 /im L = 1.3 /im and 
W = 4.8 /im, L = 0.7 /im, where L is the distance be- 
tween the voltage probes covering nearly the entire sam- 
ple width W (see Fig. 1). The contact resistances were 
measured to be 20 Q for the BL and calculated to be 200 
Q for the ML, and were subtracted. The gate voltage V g 
dependence of a is shown for both samples for a range 
of magnetic fields in Fig. 1. At zero field, one observes 
a slightly sub-linear increase of the conductance on both 
sides of its minimum at the neutrality point. The mo- 
bility varies between 3000 and 5000 cm 2 V -1 s _1 for the 
ML, 3000 and 6000 cm 2 V" 1 s" 1 for the BL. Above 2T, 
steps in the conductance of the ML occur near quantized 
values 4(n + l/2)e 2 //i, as expected. The oscillations in 
the BL (with a maximum of conductance at the neutral- 
ity point) look more unusual but can be understood given 
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FIG. 1: Gate voltage dependence of the conductance at sev- 
eral magnetic fields, 1 T apart. The contact resistances have 
been subtracted. Top panel: monolayer A. Bottom panel: 
bilayer B. Inset: electron micrographs of the samples. 



the aspect ratio of the sample (see below). 

We now describe how we extract r tr and r e from 
the magnetoresistance (MR) (see Fig. 2). The two- 
terminal MR results from mixing of the diagonal (p xx ) 
and off-diagonal (p xy ) components of the resitivity tensor 
[14j[T5] . The degree of mixing depends on the aspect ratio 
of the sample. For a square geometry, close to that of the 

1/2 

monolayer, R(B) = [p^ x + p^ y ] ; in a short wide sam- 
ple such as the bilayer R(B) = (L/W) [p 2 xx + p xy \ / ' p xx . 
Intermediate geometries can be calculated following the 
model developed in [14] . It is then possible to reconstruct 
the complete MR from the expressions of the resistivity 
tensor [16] valid in the limit of moderate magnetic field 
where ShdH oscillations can be approximated by their 
first harmonics: 



5pxx(B)/p = 4D T exp |-^r| cos 

Pxy(B) = pQUJ c T tr - S p xx (B) /2uJ c T tr , 



Hide 



(1) 



where po = l/ a * s the zero-field resistivity and uo c = 
eB/m* is the cyclotron frequency, m* = fikp/vp is the 
cyclotron mass which depends explicitly on the Fermi 
wave vector kp for the ML (constant Fermi velocity). 
On the other hand, the bilayer's dispersion relation is 
parabolic at low energy and m* can be approximated by 



the effective mass ra e // = 0.035m e , nearly independent 
of the carrier density in the range of V g explored where 
\Ep\ < 80 meV is smaller by a factor 5 than the energy 
band splitting [1 j. This value of m e ff is confirmed by the 
analysis of the temperature dependence of ShdH oscilla- 
tions [17 . The phase (/>, either tt or 27r, and the parameter 
j, either 1 or 2, depend on the nature of the sample (ML 
or BL). The Fermi energy Ep is TikpVp for the mono- 
layer and % 2 k 2 F /(2m e ff) for the bilayer. The prefactor 
Dp = 7/ sinh(7) with 7 = 27r 2 /c#T/fkj c describes the 
thermal damping of the oscillations. 

To analyze the data we first deduce kp from the pe- 
riodicity of the ShdH oscillations function of 1/B. We 
think that this determination is more reliable close to 
the neutrality point where the sample is possibly inho- 
mogeneous than the estimation of n c = k F /ir from the 
gate voltage and the capacitance between the doped sili- 
con substrate and the graphene sample [13 . Knowing kp 
we then determine r tr from the low field quadratic low 
field magnetoresistance which is found to be independent 
of temperature between 1 and 4 K: 



R(B) - R(0) 
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a g (u c r tr ) 2 



(2) 



We have used the relation a = p^ 1 = (2e 2 /h)kpVpT tr . 
The dimensionless coefficient a g , which depends on the 
aspect ratio of the sample is determined numerically fol- 
lowing [14] and the experimental values of Wand L. It 
is found equal to 0.53 ± 0.01 and 0.84 ± 0.02 for our ML 
and BL samples A and B, respectively. It is important to 
note that this determination of r tr is independent of any 
assumption of the contact resistance on our two termi- 
nal samples. We finally extract r e from the damping of 
the first harmonic of ShdH oscillations in the resistivity 
tensor in exp(— f3/B) where f3 = nhkp / 'evpr e , see Eq{l] 
The kp dependences of r tr and r tr /r e are shown in 
Fig. 3 for samples A and B as well as three other ML sam- 
ples, consisting of another two-terminal sample C (very 
similar to A) and two multi-terminal samples (D and E) 
with Hall-bar geometry (see [25] for more details). We 
observe different behaviors for the ML samples, where 
r tr has a minimum at the CNP, and the BL, where it has 
a maximum. In all cases, despite rather large variations 
of Tt r , T tr /T e is nearly independent of kp. It is equal 
to 1.7 ± 0.3 for the monolayers A,C,E and to 1.8 ± 0.2 
for the bilayer in the whole range explored, which corre- 
sponds to n c between 1.5 x 10 11 and 5 x 10 12 cm -2 . That 
r tr /T e is of the order but smaller than 2 indicates that the 
typical size of the scatterers does not exceed the Fermi 
wavelength. We note however that sample D exhibits a 
value of r tr /r e at high electron doping which is larger 
than 2 (~ 2.4). The area of this sample (12pm 2 ) is much 
larger than the area (~ lpm 2 ) of all the other samples 
A, B, C and E. We suspect that this large sample con- 
tains more spatial inhomogenities than the other smaller 
samples which could explain a reduced value of r e . 
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FIG. 2: Analysis of the magnetoresistance. Left panel: Mag- 
netoresistance of monolayer sample A. Dots: experimental 
points at T — IK; Continuous line: fit according to Eqs. |2]) 
and Inset: B 2 dependence of the low- field magnetore- 
sistance for different gate voltages (Curves shifted along the 
Y axis for clarity). rt r is extracted from the slopes of these 
curves according to Eq. |2]). Notice that the slope increases 
in the vicinity of the Dirac point reflecting the divergence of 
the inverse effective mass. Right panel: ShdH oscillations 
of the longitudinal component of the resistivity in bilayer 
sample B for different temperatures after subtraction of the 
quadratic backround. The Fermi wave vector &f and the elas- 
tic time r e are deduced from the period and the decay of the 
oscillations with 1/B at low temperature. Inset: Tempera- 
ture dependence of the oscillations amplitude normalized to 
T = 0. Solid line: fit according to the Lifshitz-Kosevich for- 
mula Dt = 7/ sinh(7) with 7 = 2n 2 ksT /Tiuj c ,IS]. The effec- 
tive mass determined from this fit is m e ff = 0.035 =b 0.002ra e 
in the whole range of gate voltage investigated. 



Finally, it is also possible to fit cr(V g ) at 5T depicted 
in Fig.l using the value of r e determined as described 
above. Taking into account the geometry of samples A 
and B, following [14], one can relate the contributions of 
the n th Landau level to the conductivity tensor, S n a xx 
and 5 n <j xy within the semi-circular model. 

5ncr x x ~ exp -ln2 [{v - (y n + v n+1 ) /2) /T v ) 2 ^ 

Here, v n is the filling factor of the n th level, v is the filling 
factor in between the n th and (n + l) th levels, and cr® y n is 
the quantized Hall conductivity at the n th plateau equal 
to 4(n + l/2)e 2 /h and 4ne 2 /h for the ML and BL, re- 
spectively. The width of Landau levels function of filling 
factor, IV, and energy Te = h^/2ujr-/7rr^fi8\. are related 
via: 

r -<n/*f- (4) 

Good agreement between the experimental data and the 
two-terminal G(V g ) of the ML, calculated from the con- 
ductivity tensor [14] [15], is obtained taking the filling 



FIG. 3: Jcf dependence of r tr and Tt r /r e ratio. Left panel: 
monolayers A, C, D and E. Right panel: bilayer B. The con- 
tinuous lines are the fits for samples A, B and D according 
to the resonant impurity model, Eq. [5] For samples A B and 
C (two terminal configuration) rt r was extracted from the 
low field magnetoresistance (crosses) whereas it was extracted 
from the zero field conductivity for samples D and E. Posi- 
tive/negative values of Jcf correspond to electron/hole doping. 
Lower panels: ratio rt r /r e where r e is deduced from the fit of 
the low temperature decay of the ShdH oscillations. Dotted 
lines figure the average value rt r /r e = 1.8. Interestingly, al- 
though the mobilities and accordingly r tr vary substantially 
from one sample to the other (from 5000 to 800 cm 2 /Vs from 
samples A to E) the ratio rt r /r e is similar for all samples. 



factor dependence of T v from Eq. Q, see Fig. 4. For 
the BL, we had to modify the semi-circular relation in a 
similar way as done in Ref. [T5] . 

We now compare our results on r e and r tr to recent 
theoretical predictions. We first consider scattering on 
charged impurities [8j[9]. In particular, the question of 
the difference between r e and r tr has been addressed for 
a graphene monolayer [19]. The minimum value of the 
ratio T tr /T e is obtained when the impurities are located 
close to the graphene foil, in which case it is expected 
to be 2 - as a result of the absence of backscattering 
- and independent of n c . This is a bit larger than the 
measured ratio for most samples. Screened charged im- 
purities are characterized by a screening radius l/q sc , 
which in the Thomas-Fermi approximation, is given by 
1 / Qt f = 7rehvF/e 2 kF, where e is the appropriate dielec- 
tric constant. In the Born approximation, the transport 
time is r tr ex c^ sc vfI^f- For a monolayer, ^tf/^f is a 
constant ~ 3 and both r tr and r e are then expected to 
increase as hp, which is not what we observe in Fig. 3 
where the increase is sublinear. The disagreement is even 
stronger for a bilayer, where the ratio ^tf/^f oc 1/kp 
varies between 3 at high doping and 12 close to the neu- 
trality point. The transport time is then expected to 
vary linearly with n c , if the screening radius is estimated 
as ~ 1/&f, or to be independent of kp if estimated as 
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~ 1/qrF <^ 0, neither of which agrees with our 

data, see Fig. 3. 

An alternative explanation is resonant scattering re- 
sulting from vacancies or any other kind of impurities of 
range R such that a ^ C where a is the carbon- 

carbon distance, and with a large potential energy [7tfT2]. 
It is characterized by a transport cross section 



k F \n 2 (k F R) 



(5) 



The resulting transport time r tr = l/(niV F A tr ) (n^ is 
the concentration of impurities) leads to a conductance 
increasing as n c with logarithmic corrections for both the 
ML and BL. In both cases, our extracted r tr (kF) (see 
Fig. 3) are compatible with the square logarithmic de- 
pendence of Eq. ([5]). It is also possible to estimate the 
range of the impurity potential 0.5 A < R < 2. 5 A and 
the concentration of impurities rii = (8 ± 2)xlO n cm -2 , 
which turns out to be identical for sample A and B. This 
is of the order of the minimum value of the carrier den- 
sity rimin = 1.5 x 10 11 cm -2 , extracted from the exper- 
iment. It is also interesting to note that the minimum 
conductivity expected for this resonant impurity model, 
crmin = (2e 2 /nh)(n min /ni)ln 2 (Ry/nn min ) = 3.7e 2 //iand 
4.5e 2 /h for the ML and the BL, respectively, are similar 
to the observed experimental values which are 3.3 and 
4.1 e 2 /h. This analysis also corroborates our results on 
the ratio r tr /T e indicating scatterers with a range smaller 
than the Fermi wavelength (but possibly of the order if or 
slightly larger than the lattice spacing) . Whereas the res- 
onant character is not essential for the validity of Eq. ([5| 
for massive carriers (corresponding to the bilayer) [20] , it 
has been shown that it is essential for massless carriers 
in the monolayer [21 . This resonant-like character, al- 
though not straightforward, has been demonstrated in 
the case of scattering centers created by vacancies in 
graphene over a wide range of Fermi energies [22 . As 
shown in detail in [25 , it is not necessary to fine-tune kp 
to obtain the In 2 dependence in Eq. |H]). 

In conclusion, our results indicate that the main scat- 
tering mechanism in our graphene samples is due to 
strong neutral defects, with a range shorter than the 
Fermi wavelength and possibly of the order of a, induc- 
ing resonant (but not unitary) scattering. Likely candi- 
dates are vacancies, as observed recently in transmission 
electron microscopy [23] , voids, ad-atoms or short range 
ripples as suggested in [23] . This does not exclude the 
presence of long range charged impurities responsible for 
electron hole puddles but their contribution to the scat- 
tering rates l/r tr and l/r e appears to be negligible in all 
the samples investigated. 
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FIG. 4: Comparison of G(v) at 5T for samples A and B with 
the expression of the conductance derived in Ref. [14], taking 
the aspect ratio L/W of each sample and using Eqs. |3| and 
Q with r e (/ci?) determined above. The dashed vertical lines 
indicate the positions of v n which, as expected, are different 
for the mono v n — ±4(n +1/2) and the bilayer v n — ±4n . 
The conductance quantization is well obeyed for the mono- 
layer but not for the bilayer. This is well explained by the 
aspect ratio of the bilayer sample [21 US] . 
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